Superfluid-insulator transition in a periodically driven optical lattice 
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We demonstrate that the transition from a superfluid to a Mott insulator in the Bose-Hubbard 
model can be induced by an oscillating force through an effective renormalization of the tunneling 
matrix element. The mechanism involves adiabatic following of Floquet states, and can be tested 
experimentally with Bose-Einstein condensates in periodically driven optical lattices. Its extension 
from small to very large systems yields nontrivial information on the condensate dynamics. 
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The Bose-Hubbard model plays an important role in 
condensed matter physics, since it embodies essential fea- 
tures of strongly interacting Bose systems in a minimal 
manner, namely the competition between kinetic and po- 
tential energy effects, and the resulting quantum phase 
transition from a superfluid to a Mott insulator Q, 
It describes Bose particles on a lattice with on-site in- 
teraction, so that particles occupying the same lattice 
site repel each other, while tunneling is allowed between 
adjacent sites. This is expressed by the Hamiltonian 
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where cp is an annihilation (creation) operator for a 

boson on the site labeled j, and hj — CjCj denotes the 
corresponding number operator. The first sum runs over 
all pairs of neighboring sites i and j, with the matrix ele- 
ment J quantifying the strength of the tunneling contact. 
Moreover, U is the repulsion energy contributed by one 
pair of bosons located on the same site. Therefore, the 
characteristic dimensionless parameter is the ratio U/J: 
When U <C J, so that tunneling dominates, the ground 
state of the system describes a superfluid, whereas it has 
the properties of a Mott insulator when the interaction 
dominates, U ^> J . Mean-field theory 0,0 gives the crit- 
ical value (U/J) c «zx 5.83 for the transition in a lattice 
filled with one particle per site, which captures the case 
of a three-dimensional (3d) cubic lattice with coordina- 
tion number z = 6 reasonably well, whereas more refined 
methods yield a Kosterlitz-Thouless transition with 
(U/J) c « 3.8 for d= 1. 

After the model QJ had long been of primarily theoret- 
ical interest, it has found its laboratory realization with 
Bose-Einstein condensates in optical lattices jj, |5j, |6J . In 
such systems, the expected transition has been observed 
upon varying the lattice depth, both for d = 3 and 
d = 1 8j. In this Letter, we demonstrate that the transi- 
tion from a superfluid to a Mott insulator can be induced 
in an altogether different manner which, in contrast to all 
scenarios studied before, hinges on the effect of a time- 
dependent force, and which can be assessed experimen- 
tally with condensates in periodically modulated optical 



lattices. We will investigate a periodically forced Bose- 
Hubbard model for d — 1, as described by the explicitly 
time-dependent Hamiltonian 

H{t) = H +K cos(ujt) ^2 j fij , (2) 
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where the equidistant sites are labeled according to their 
position in ascending order. The oscillating term, which 
mimics a monochromatic electric dipole potential with 
frequency u> and amplitude K, can be realized experi- 
mentally by periodically shifting the position of a mir- 
ror employed to generate the standing laser wave, and 
transforming to the co- moving frame of reference |y, [lOj . 
We will argue that the driven system behaves, for 
sufficiently high frequencies, similar as the undriven sys- 
tem Q , but with the tunneling matrix element J of the 
latter being replaced by the effective matrix element 

J eff = J3 {K/(huj) , (3) 

where Jo(a^) denotes the ordinary Bessel function of or- 
der zero. Hence, the actual control parameter becomes 
U I J c g, which can be varied by adjusting the parameters 
of the periodic modulation. This implies the possibility 
to switch between the superfluid and the insulator state 
by changing, e.g., the modulation strength K . 

The rescaling (J2J) is not unfamiliar with periodically 
driven single-particle quantum systems. It occurs, among 
others, when a particle moves on a periodically forced 
Id lattice with nearest neighbor coupling 0, Hp. such 
as an electron in a semiconductor superlattice |l3j |. It 
also underlies the Jo-type renormalization of atomic g- 
factors in oscillating magnetic fields 0, 0], and the 
coherent destruction of tunn eling of a particle in a pe- 
riodically forced double- well [161 1 1 7| . However, as will 
be discussed below, the many-body system |j5J| is signif- 
icantly more involved when the thermodynamic limit is 
taken; the rescaling © then describes only part of the 
relevant physics. 

Our analysis is based on quantum Floquet the- 
ory ^^|: Since the Hamiltonian © depends periodi- 
cally on time, H(t) = Hit + T) with period T = 2?r/w, 
there exists a complete set of solutions to the time- 
dependent many-body Schrodingcr equation of the form 
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|^n(i)) = \u n (t)) exp(— ie n t/h), where the Floquet func- 
tions \u n (t)} inherit the period of the driving force, sat- 
isfying \u n (t)) = \u n (t + T)). Thus, Floquet states for 
periodically time-dependent quantum systems, obtained 
by solving the eigenvalue equation 



(H(t) - ihdt) \u n (t)) = e n \u n {t)) , 



(4) 



constitute an analog of Bloch states known from spa- 
tially periodic crystals; the eigenvalues e„, which describe 
the time evolution of these states in close analogy to 
the evolution of energy eigenstates, are called quasiencr- 
gies. While in solid-state physics quasimomenta are de- 
fined up to an integer multiple of a reciprocal lattice 
vector, quasienergies are defined up to an integer mul- 
tiple of hw: If \u n {t)) solves Eq. (@J with eigenvalue s n , 
and m = 0, ±1,±2, then \u n (t)) exp(imwi) is a T- 
periodic eigensolution with quasienergy e n + mhuj. The 
quasienergy spectrum of a periodically time-dependent 
quantum system thus possesses a Brillouin zone-like 
structure, the width of one zone being hw. 
We then employ the Floquet basis 



\{nj},m) = |{n j })exp 



l- — sm(ti)t) y jrij 



hw 
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lmujt 



(5) 

where indicates a Fock state with rij particles on 

the jth site, and m again accounts for the zone structure. 
The eigenvalue problem refers to an extended Hilbert 
space of T-periodic functions, in which the time variable 
is regarded as a coordinate [ijj , so that the scalar product 
in that space is given by 
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dfct 



(6) 



i.e., by the usual scalar product (-|-) combined with time- 
averaging. Hence, the quasienergies are obtained by com- 
puting the matrix elements of the operator H (t) — ihdt 
in the basis J5j with respect to the scalar product JBJ, 
and diagonalizing. Denoting the Hamiltonian Hq with 
J = 0, which is diagonal in the basis by H ln t, and its 
J-proportional tunneling term by H tU m we find 



= S„ 



({nj}\Hi n t\{nj}) + mhw 



l J m ^ m (K/(huj)) ({n^l-HtuaKni}) , (7) 



where s = J2j( n 'j ~ n j)j = il> smce -Htun only transfers 
one particle by one site. We observe that with respect to 
the "photon" index m this matrix has a transparent block 



structure: The diagonal blocks with m 



reproduce 



the matrix which yields the eigenvalues of the undriven 
system, but with J replaced by J Q g according to Eq. ||3J| , 
and replicas shifted by integer multiples of hw. These 
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FIG. 1: Exact energy spectrum of a small one-dimensional 
undriven lattice £Q| with M = 5 sites and N = 5 particles, 
versus Mott-Hubbard parameter U/J. The splitting-off of 
the ground state with increasing interaction strength can be 
regarded as a precursor of the quantum phase transition. 



blocks are coupled by nondiagonal ones proportional to 
Bessel functions Z m i^ m yK/{hS)j. Obviously, the antic- 
ipated rescaling 10 holds only to the extent that these 
couplings can be neglected. This will be the case, at least 
in a perturbative sense, if the block separation hw is much 
larger than both the energy scale J of the coupling and 
the energy scale U associated with the diagonal blocks, 
i.e., for high frequencies hw 3> max{ J, U}. 

To demonstrate that this reasoning is justified, we 
present numerical results for small systems. Figure Q 
depicts the exact energy spectrum for a one-dimensional 
undriven model with N = 5 particles on M = 5 
sites. Even here, the precursor of the superfluid-insulator 
transition already is apparent: With increasing U/J the 
system's ground state, associated with a uniform distri- 
bution of the particles over the sites, splits off from the 
group of excited states, which describe various patterns 
of particle-hole excitations. In the limit of an infinitely 
large system, N — > oo and M — > oo with N/M = 1 
held constant, the excited states form continuous energy 
bands; the ground state then splits off from the lowest 
band at a finite (U/J) c |3|]. This separation of an indi- 
vidual state from the continuum indicates the transition 
to the Mott insulator state. 

For comparison, Fig.[5]shows the first Brillouin zone of 
numerically computed quasienergies for the driven sys- 
tem (J2J, again with N = M = 5, scaled frequency 
hw/ J = 14, and scaled driving amplitude K/hw = 1.5. 
Since Jq(1.5) k 0.5, this set of parameters allows for a 
convenient test of the hypothesis ©: With J e g « J/2, 
the quasienergy spectrum of the driven system J3J for 
a given parameter U/J should correspond (apart from 
its zone structure) to the energy spectrum of the un- 
driven system Q with the same U/J c g, which is about 
2U/J. This is borne out, to remarkable accuracy by 



FIG. 2: First Brillouin zone of the exact quasienergy spectrum 
of the driven Id system with N = M = 5, huj/J = 14, 
and K/huj = 1.5. This spectrum is recovered approximately 
from the energy eigenvalues shown in Fig. if J is replaced 
by Jeff « J/2 here, and the energies are taken modulo hu). 

a comparison of Figs. ^ and |2 With the scales of the 
respective axes differing by a factor of two, the eigenval- 
ues plotted in Fig. 2 almost equal those in Fig. 1. As a 
consequence of the Brillouin zone structure, quasienergy 
eigenvalues which disappear at the upper zone bound- 
ary reappear again at the lower one. This reappearance 
is a source of substantial complications: States originat- 
ing from different Brillouin zones are coupled through 
the matrix elements neglected in the explanation of the 
renormalization JSJ), so that many apparent level cross- 
ings in Fig. factually are tiny avoided crossings. It is the 
high-frequency condition fko 3> max{J, U} which guar- 
antees that these avoided crossings remain too narrow 
to be resolved. If this condition is not met, a multitude 
of large avoided crossings appears in the spectrum, thus 
revealing typical signatures of quantum chaos |2ij • 

Figure |3 shows quasienergies for U/J — 3 kept fixed, 
again for Hlj/ J — 14, as functions of the scaled amplitude 
K/(j\S). For K/{fajS) « 2.4, close to the first zero of Jo, 
the tunneling contact is quenched almost entirely, so that 
the various bands of particle-hole excitations collapse. In 
the vicinity of this point, the Floquet state evolving from 
the unperturbed ground state must have the properties 
of a Mott insulator, although the ground state of the 
undriven system describes a superfluid for U/J = 3. 

To verify this conclusion, we have solved the time- 
dependent Schrodinger equation for a system with N — 
M = 7 and U/J = 3, initially prepared in its superfluid 
ground state and then subjected to periodic forcing with 
frequency fko/J — 14, and an amplitude K/(hio) which 
increases linearly from to 2.4 during the first 100 cy- 
cles T, then stays constant for another 100 cycles, and 
finally is linearly ramped down to zero between t = 200 T 
and t = 300 T. Experimentally, the superfluid phase is 
detected by a sharply peaked reciprocal lattice pattern in 



FIG. 3: First Brillouin zone of the exact quasienergy spec- 
trum of the driven Id system with N = M = 5, hu>/ J — 14, 
and U/J = 3, versus scaled driving amplitude K/hui. For 
K/hu) ~ 2.4, close to the first zero of Jo, the tunneling con- 
tact is (almost) switched off, resulting in a collapse of the 
different bands. In the vicinity of this value, the Floquet 
state evolving from the system's ground state (marked by the 
arrow) has the properties of a Mott insulator, even though 
the undriven system's ground state is superfluid for U / J = 3. 
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FIG. 4: Time evolution of the momentum distribution g(p) 
obtained by solving the time-dependent Schrodinger equation 
for N = M = 7, hcu/J = 14, and a slowly varying amplitude 
Kj (huj) which is linearly ramped up to 2.4 during the first 100 
cycles T, then stays constant for another 100 T, and is ramped 
down to zero during the final 100 cycles. The system was 
initially in its ground state; the distribution was recorded at 
integer multiples of T. The disappearance and reappearance 
of the peak pattern signals the transition from the superfluid 
to the Mott state and back. 



the momentum distribution, g(p) = J^j fc( c j c fc) e_lp ^~ fc \ 
which can be measured by time-of-flight absorption imag- 
ing. (Here, p is given in multiples of ft/lattice constant.) 
Figure ^ depicts our result: Since Floq uet states respect 
an approximate adiabatic principle [2l| , the initial super- 
fluid ground state is first adiabatically transformed into a 
Mott insulator state, as witnessed by the disappearance 
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of the peaked momentum distribution, and then trans- 
formed back to the inital state, apart from a remaining 
excitation of other states totalling to a few percent. This 
is a major result: The amplitude of the periodic force 
decides whether the system is superfluid, or in a Mott 
insulator state. 

When extrapolating from these model calculations to 
large systems, two issues have to be considered. Firstly, 
when increasing the number of lattices sites while main- 
taining an occupancy of one particle per site, say, the 
quasienergy levels fill the Brillouin zone densely; in the 
thermodynamic limit, the spectrum probably is a contin- 
uum. Then there will be no "sharp" Floquet state evolv- 
ing from the ground state, but rather a resonance with 
a finite lifetime, due to the residual couplings to other 
states. Starting from an undriven, infinite system with 
a superfluid ground state, and switching on the periodic 
force, we conjecture that a Mott-msulator-like resonance 
appears at that amplitude K which, after rescaling ac- 
cording to Eq. , corresponds to that tunneling matrix 
element J which marks the quantum phase transition in 
the undriven system. Secondly, for an infinitely large 
lattice there is no adiabatic limit when switching on the 
driving force [12]; turn-on and turn-off necessarily have 
to take place within a short interval. Our calculations in- 
dicate that adiabatic following can even be improved by 
shortening the turn-on time, since Landau-Zener transi- 
tions at narrow avoided crossings, possibly corresponding 
to condensate heating, then are suppressed. 

While these issues are still not covered by rigorous 
mathematical theorems on periodically driven quantum 
systems |23j , and remain out of reach of even most pow- 
erful supercomputers, they can be addressed in the lab- 
oratory. In experiments with cold atoms in driven op- 
tical lattices [l(j|, narrowing of Bloch bands compatible 
with the rescaling has already been observed, even 
though the single-band regime has not been reached. 
Employing Bose-Einstein condensates in optical lattices 
in order to realize the driven Bose-Hubbard model J5J, 
one has to respect not only the high-frequency condition 
Huj > max{ J, U} required for the approximate Bessel- 
function rescaling ©, but there is the obvious additional 
condition ftw < A, where A denotes the gap between the 
lowest two Bloch bands of the undriven lattice, in order 
to exclude transitions to higher band states. Elementary 
estimates in the spirit of Ref. [5j suggest that under typ- 
ical conditions (as provided by 87 Rb atoms in a lattice 
created by laser radiation of A ~ 852 nm wavelength 0) 
this leaves a viable window of frequencies in the low 
kHz regime. For higher filling factors, or in 3d lattices 
with forcing in all three directions, the critical parame- 
ter (U / J) c becomes much larger, allowing one to employ 
deeper lattices with larger band gap A, and hence to 
work with still higher frequencies without violating the 
single-band approximation. A quantity of key interest in 
such experiments will be the extent to which, after start- 



ing from a superfluid ground state, then ramping up the 
force into the insulator regime and ramping it down again 
as in Fig. 01 the superfluid peak pattern reappears, pro- 
viding information on both the lifetime of the conjectured 
Mott-like resonance state and the degree of adiabatic fol- 
lowing, or, more generally, on the extent to which the 
quantum evolution of a mesoscopic matter wave can be 
guided even under critical conditions. Thus, the scenario 
envisioned here is not intended as a look at the common 
superfluid-insulator transition from a different angle, but 
aims at obtaining genuinely new, nontrivial information 
on condensate dynamics. 
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